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I. Introduction 



The events of ultra-high energy cosmic rays (UHECR) corresponding to 
primary energy above 10 19 eV are difficult to explain© with conventional as- 
trophysical objets, both regarding their acceleration mechanism and propa- 
gation towards the earth due to the interaction with the cosmic background 
radiation! (CBR) if the source is beyond ~ 50Mpc. 

A possible solution of this enigma is given by the so called top-down 
mechanism! where long-lived very massive microscopical objects decay pro- 
ducing the UHECR, which is plausible because so far the events of the latter 
appear to be roughly isotropic. 

In any case the top-down mechanism would imply physics beyond the 
standard model of particles and their interactions (SM). One alternative cor- 
responds to superheavy relics! , quasi-stable because their interactions with 
the known particles are of gravitational order, which might belong to the 
hidden sector where supersymmetry is broken. Another possibility is given 
by cosmic strings^ formed in the phase transition due to the breaking of a 
symmetry at the scale of a grand unification theory (GUT). Though the or- 
dinary Kibble strings! consisting of Higgs and gauge fields might explain! 
the UHECR using their flux to normalize the model, either the insertion of 
monopoles forming necklaces^ or the superconducting Witten strings! with 
the addition of fermionic fields which give quasi-stable closed loops called 
vortons seem more suitable. 

It is the purpose of the present work to analyze the details of vortons 
as source of UHECR. It had already been seen^ that they may produce 
the global flux above 10 19 ey provided their density is dramatically reduced 
during the electroweak transition in a way that might generate the matter- 
antimatter asymmetry of the universe0 . 

We now describe the energy spectrum of the flux given by the decay 
of a superheavy boson emitted by the vorton and conjecture, because of 
the hadronization of the resulting quark, that it can constitute the hard 
component that might emerge above the GZK cutoff according to the recent 
presentation of events^! . Even though vortons should behave as cold dark 
matter (CDM) and be concentrated in the galactic halo, we also study the 
softening caused by redshift if they were distributed uniformly in the space as 
occurs for necklaces and find that it is not possible to distinguish both cases 
at present, though in the latter situation a more important effect coming 
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from the interaction with CBR should produce a depression due to the GZK 
cutoff followed by a recover of the spectrumEl caused by its hard nature at 
emission. 

We then face the difficult problem of the dynamics of strings, crucial for 
determining the vorton density before and after the electroweak transition. 
We assume that the GUT model allows superconductivity to appear at the 
same scale of the string formation and evaluate that the delay in the sta- 
bilization of vortons is not too relevant to reduce their density. When the 
universe cools down to the electroweak temperature, we estimate the condi- 
tions for the rate of destabilization of vortons due to disappearance of the 
Dirac zero-modes which originated the superconductivity of exotic quarks 
and the replacement with those of ordinary fermions. The result is that the 
collapse of the most abundant short vortons is sufficiently gradual to avoid 
the reheating that would dilute the baryogenesis, and the surviving long ones 
succeed in absorbing the ordinary fermions with parameters compatible with 
the density and lifetime necessary to explain the UHECR. 

Finally we discuss which is the possible GUT group consistent with our 
mechanism. We see that SO(10) is not adequate since at GUT scale only 
vortons with ur might be formed, which could not give the baryogenesis 
at the electroweak transition where it is moreover unlikely that long loops 
might be stabilized by ordinary fermions. On the contrary, E 6 is suitable for 
our purposes because at high temperature vortons with exotic quarks may be 
formed linked to zero-modes that subsequently disappear at electroweak scale 
due to mass effects of the light Higgs, whereas new zero-modes for ordinary 
fermions are allowed by the existence of two additional abelian symmetries 
of the model apart from the electromagnetic one. 



II. Energy spectrum of UHECR flux from vorton de- 
cay 

Considering sources that emit n (£) UHECR per unit space and time the 
total flux on earth will be 
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where a is the scale parameter of universe, to its present age and ti n an initial 
time which depends on the assumed distribution of sources but in no case is 
earlier than that which by redshift produces energies at least ~ 10 19 eU. 

If the sources are quasi-stable objects like vortons with density n(t) each 
one having at time t a probability T per unit time of emitting UHECR 
observed on earth, 

F= ±j dtn{t) (m 3 r . ( 2) 

We assume that the vorton emits by tunneling a superheavy particle X, 
Higgs or gauge boson of GUT scale, which very quickly decays in quarks and 
leptons, the former giving the UHECR by hadronization. We will consider 
two cases: the more plausible one in which vortons, behaving as nonrelativis- 
tic particles, condense in the galactic halo, and the other extreme alternative 
analogous to ordinary cosmic strings in which they are still uniformly dis- 
tributed in space. 

A. Condensation in halo 

In this case redshift may be neglected and T = — , where r is the vorton 
lifetime for emission of X whose decay produces N c UHECR. Then the total 
flux will be 

N c , x At 

F h = -^n h (to)— , (3) 

where At ~ 50kpc due to the halo size. As it will be seen in the next section, 
with r larger than t and n^(to) a fraction of the dark matter in the galactic 
halo, a flux of UHECR of the expected order is obtained. 

We now turn to the energy spectrum F(E) such that F — / dE F(E), 
where the limits of integration correspond to the UHECR range. From the 
probability distribution 

dT 1 Nc 

TE = -^(E-E t) , (4) 

the flux spectrum will be 

1 At Nc 

F h (E) = n h (t )—J2S(E-E t ) . (5) 
4tt r i=1 
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To compare with observations, one must average on the intervals AEi 
separating neighbouring particles in energy 

Fh{El) ~ A~E~ "Itt ~r ' (6) 

If events are equally spaced in log E, which is plausible in hadronization 
with QCD except for the upper limit ~ m x as it will be discussed below, 

AEi ~ Ei , ~F h (Ei) ~ -1 , (7) 

Ei 

corresponding to a hard component compared with the standard behaviour 
for lower energy which is roughly F(E) ~ E~ 3 . 



B. Uniform distribution in universe 

Since we consider vortons as quasi-stable particles, if we assume a hypo- 
thetical uniform density n u (t Q ) at present, from Eq.(2) the spectrum at earth 
would be 

F^^L^iy^-^, ( 8) 

where the lower limit is approximatively the matter-radiation equivalence 
time t eq in order to include into UHECR particles redshifted from ~ 10 2i eV 
which we take as the maximum energy. If we disregard for the moment the 
attenuation due to interaction with CBR, the probability distribution is the 
same as Eq. (4) but referred to emission energy^ E em 

dT 1 Ne 

= -J2S(E em -E t ) , (9) 



dE em , ■ _ i 



where the relation with observed energy is given by redshift z 



t \ 9 



E em = (l + z)E , l + z= (^) . (10) 



Therefore the spectrum Eq.(8) becomes, being t » t 



eq i 



Q + N c -1 

F U (E) = —n u (t )^E- 2 J2 —OiE, - E) . (11) 
87r r i=i Ef 
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The total flux for UHECR denned for E > E is, with E < E { 



roc | f n 

F u = dE F U {E) = — n u (t )- 

JEo 47T r 



Averaging the spectrum in intervals defined for convenience as 



F u (Ej 



AEj Je^ 

and with the above hypothesis of AEj ~ Ej 



FJEi) = — njt ) — — — 



dE F U {E) 



E 



;i2) 



(13) 



(14) 



it is clear that redshift will produce a softening of the law 4- through the 
bracket factor, which will be slight except for the highest j. 



C. Comparison with observations 

Considering that the probability per unit time for a vorton to emit par- 
ticles with energy between E^ — 10 19 ey and Eh — 10 2A eV from the decay 
of X is 

r E H dT N 

(15) 



> El dE em 



dEf,„ 



Nr. 



T 



the definition of the average on intervals according to Eqs. (9) and (7) 

dT _ 1 1 
dE ~~ ~r~E ' 



(16) 



gives, from Eq.(15), N c ~ 10 which is a reasonable value compared with the 
extrapolation of fragmentation function^ . This is consistent with the total 
energy emitted by a vorton per unit time, being nix ~ E H , 



e h dT 1 

dE em —— — E em — —E H 
o dE pm t 



(17) 



Therefore we may take the equally spaced particles in log-E according 
to Ei ~ 10 19 ey, E 2 ~ 10 19 - 5 eV . . . E 9 = 10 23 eV, E 10 ~ 10 23 5 eV, and 
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E ~ 10 18 - 5 eV, so that with the previous definition AEj = (l — y^ts) Ej 
quite compatible with Eq.(16). 

For the vortons condensed in the halo, since from Eq. (6) the flux in each 
bin is the same, one might roughly normalize it at the observed value for 
10 19 eV i.e. 

n h {t ) At = 1 
47T r km 2 yr 

Being the mass of a vorton ~ N^mx and representing a fraction / of the 
halo average energy density ~ 0.3 it turns out 

n h {t ) At ^ f t 10 7 
47r t Nl t km 2 yr ' 

so that for r ~ to and Nl ~ 10 3 which is a sensible number as will be 
discussed in Sec. Ill , a fraction / ~ 10~ 4 of dark matter would be enough. 
More precisely, from the recent presentation of data shown in Fig.l, there 
seems to be an extragalactic component above the "ankle" and then the hard 
component dominating beyond the GZK cutoff. With this interpretation the 
contribution of vortons is a small part of flux at 10 19 eV and / may be even 
two orders of magnitude smaller. The fact that our fit normalizes the hard 
component at ~ 10 20 eV allows to reproduce data with nix ~ 10 15 GeV at 
variance with a similar discussion for superheavy relics^ . 

If we instead imagine vortons uniformly distributed, they should consti- 
tute a fraction of the critical density of universe p c (t ) ~ lCT 29 ^-. Disre- 
garding the redshift depression and GZK cutoff J = ^n u (to)^ gives a flux 
3 times larger than that of Eq.(19) because the smaller n u is compensated 
by the larger t Q compared to At. The decrease due to redshift corresponds 
to the bracket in Eq.(12) which is just 0.98. Similarly the bracket in the 
energy spectrum of Eq.(14) gives F U (E) oc E( i +K ) where K ~ 0.011 as is 
seen by the fit of Fig. 2. Therefore it would not be possible to distinguish 
the two cases with the present statistics by the redshift effect. But certainly 
more important is the effect of the GZK cutoff which is not so drastic for the 
hard vorton componentEl giving an effective spectrum F U (E) oc 4^ between 
10 19 and 10 20 ey and a recovering behaviour oc ^3 i n the range from 10 20 
to 10 21 eV! Therefore the fraction / of dark matter should be one order of 
magnitude larger than in the case of condensation in halo. Similarly, it has 
been evaluatedliZlthat the hard component of UHECR emitted from observed 
galaxies avoids the GZK cutoff. 
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Regarding the absolute contribution of other galaxies apart from ours, 
since luminous matter is Ql ~ 5 x 10~ 3 and halos are 10 times larger £\ ~ 5 
x 10 -2 , a fraction of vortons ~ 10~ 6 would give 5 x 10~ 8 of critical density. 
This would be 200 times smaller than the required extragalactic vortons as 
said above. This estimation of the negligible contribution of halos of other 
galaxies is consistent with detailed computations!!! . 

More statistics is needed to test the anisotropy in favour of larger mass 
concentration in the halo case. It is interesting that UHECR above 5xl0 19 ey 
seems to showS a small anisotropy not related to the structure of the local 
universe. One must remark that the uniform distribution of superheavy relics 
is critically constrained by the diffuse gamma flux at GeV scale from EGRET 
requiring a low extragalactic magnetic field ~ lO~ 12 Gauss. 

A hard component similar to Eq.(7) appears from accurate calculation^ 
with QCD but it is also suggested by semiquantitative arguments. On one 
side the quark which comes from decay of X may be considered as a ultra- 
relativistic particle suffering a constant force due to friction. Therefore the 
decrease of its momentum is proportional to time and since this corresponds 
to hadronization, to produce a more energetic particle it takes more time, i.e. 
the law of Eq.(16) follows. In another way, the hadronization results from 
the emission of a gluon with energy similar to that of the quark and close to 
its direction. Thus the transition amplitude is ~ 4 and the number of final 
states EdE to keep a fixed angle around the quark, i.e. the probability of 
hadronization from vorton 

dT ~ -^-EdE , (20) 
T E 

is Eq.(16). 

All what said in this Section is valid both for vortons and superheavy 
relics. The next one will be devoted to dynamics of vortons to discuss in 
which way they may have the required density for UHECR. 



III. String dynamics and vorton densities 

The dynamics which may lead to the vorton density necessary for UHECR 
consists of several stages. First of all we will discuss the possible vorton den- 
sity above the electroweak (EW) transition temperature. For the formation 
of vortons there are four temperaturesE9 : production of ordinary strings 
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Tx, appearance of superconductivity T a , incorporation of fermionic carriers 
by loops giving the protovortons Tf, elimination of excess energy relaxing to 
classically stable configuration T r . We will take the alternative that supercon- 
ductivity, i.e. zero-modes of x-y Dirac equation in presence of bosonic string 
fields, appears in the same phase transition where strings are generated by 
Kibble mechanism i.e. T a = TV. This allows the possibility of producing the 
matter-antimatter asymmetry^ through the subsequent elimination of most 
of vortons. The number density of protovortons formed at the temperature 
T f iM 

n(T f ) ~ o , (21) 

where £ is the length below which smaller scale structure will have been 
smoothed by friction damping. Once vortons are formed, their density evo- 
lution corresponds to quasi-stable particles in an expanding universe and is 
related to that of protovortons which originate them by 




n v (T) = n(T f ) - . (22) 



If vortons could be formeclll in the friction stage of Kibble strings produced 
in a phase transition of GUT scale Tx , £ is a sort of average between the 
string damping time ~ and the Hubble time H~ x ~ i.e. 



t*yjT d H-i = {m t $%. (23) 

1 2 

Since we take T a = Tx , this will be true if the time necessary for fermions 
to be absorbed by the string forming protovortons is short. Moreover the 
subsequent interval for getting rid of excess of energy to reach the optimum 
radius of stabilized vortons must preserve the friction regime T > ^ to 
avoid radiation from the not yet static protovortons. If this occurs, from 
Eqs.(21-23) 

and since as we will show Tf is close to Tx the large density 

3 

n„(T)^pOV (25) 
\m p ij 
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is enough to produce the expected matter-antimatter asymmetry if most of 
vortons collapse at the EW transitionS . In the chiral case which corresponds 
to our fermionic carriers, their number in the loop is0 

N « iT x (26) 

and from £(T f ) of Eq.(23) N « which for the GUT scale T x ~ 

10 16 GeV will give the most abundant vortons with iV ~ 10 carriers. 

According to some approchesi! , it is hard that vortons can be formed 
in the friction stage of string dynamics. We will use a simplified model to 
estimate Tf and T r to verify that they are in the friction regime. Adopting the 
phenomenological criterium of seeing whether the fermions can be absorbed 
by the Kibble string to give way to a protovorton before the loop collapses, 
the rate of their incorporation may be given by 

~dT = an ° ' ( 27 ) 

where n D is the outher fermionic density. We assume, subject to consistency, 
that the process is fast enough to disregard the universe expansion and con- 
sider the above fixed number of fermions to be incorporated to the string. n 
will roughly correspond to a radiation mode at GUT scale n ~ T x . We re- 
quire that the density inside the string passes from zero at the time of Kibble 
string formation t = , to a final value in one direction due to field fluctu- 
ation such that niL P 7j7r ~ N , where L p is the protovorton length and ^ its 

width . The probability per unit time of fermion absorption a = aoVl ~ v 2 
will include the Lorentz factor for time dilatation due to the loop velocity 
of contraction v which is qualitatively consistent with the statement that 
vorton formation is more difficult for large velocity^ . The probability in 
the rest frame must be = h mx because it increases with the difference 
between the mass of the fermion outside due to symmetry breaking at GUT 
scale and its zero value inside the string which favours its flow there, h is a 
free parameter presumably smaller than 1 if the mass of the exotic fermion 
is smaller than that of the GUT gauge boson as occurs for most of ordinary 
fermions compared to EW bosons. Therefore we require from Eq.(27) 

N . r At 



L p 



h m x / dty/l - v 2 , (28) 
Jo 
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for At smaller than the collapse time of string . Since the protovorton must 
subsequently lose the excess of energy to reach stabilization as a vorton of 



length L v = — , the original length of the string will be L > L p > L v . 

For the step between the formation of the Kibble string of length L up 
to the absorption of N carriers to have a protovorton of length L p we consider 
that the loop contracts due to the tension [i = m 2 x , which is also the energy 
per unit length, according to 

- m\ = j t (E v) , (29) 
with E = L m x , v = it • The result is 

f = — ^ , (30) 



which, inserted into Eq.(28), gives 



£ = (31 » 



7T— 2ar sin 



where \ = , I = f dz (cos z) 2 . 

The procedure is to choose a value of h and iV = 10 as said above if 
Eq.(25) holds!. Assuming an initial length Ao, one obtains from Eq.(31) 
and then the protovorton velocity v p from Eq.(30) which must be consistent 
with the subsequent step to determine Ao . 

This corresponds to the delay between protovorton and vorton. When 
the former is planar it has an energy 

N 2 

E = ^L+—, (32) 

where, to the Kibble contribution of the first term, the kinetic one for mass- 
less fermionic carriers is added. Vortons are the classically stable loops cor- 
responding to the minimum of E. 

But on top of the energy E of the ground state of the string Eq.(32), for 
the protovorton one may add the energies corresponding to the possibility of 
twisting its pieces. The energy of each twist may be taken as 

Et = ej^j- , (33) 



11 



because for L » L v it is reasonable that it corresponds to the current, e 
being the charge. When L — > L v , E t — > oo, indicating that only the plane 
state of the vorton is possible. As a consequence, the thermodynamical 
calculation of minimization of total free energy F tot = E + F must be done 
with F = —T In Z in terms of the partition function 

N 

Z-E®^. (34) 

m=0 

where m is the number of twists. F tot (L) is a function flatter than E(L) but 
with the same minimum. Its variation is given by 

- dF tot (L) = d (Ev) v + SdT - pdL , (35) 

where ju is the chemical potential for the protovorton but, being the last term 
equivalent to — dE, the equation to be solved involves in the left hand side 
only the variation of the partition function of twists. Since moreover the 
second term on the right hand side is small because it corresponds to the re- 
heating produced by the disappearance of a piece of protovorton transformed 
into radiation, being at this stage the density of protovortons much smaller 
than that of the radiation, Eq.(35) reduces essentially to 

T dlnZ(L) ~ d[E(L)v] v . (36) 

Eq.(36) is solved numerically with the condition that at the end, when 
the stabilized vorton is reached, v — 0. Taking e ~ 0.3 thejmrtition function 
is evaluated for the two contributions N E t < T and (iV + l) E t > T as 
Z = Z\ + Z2 with 

zi =so^>^(i)(i-4 +••■)• w 

m=0 m=0 

N „ 

Z,= E = + W (37) 

m=A r +l 

Thus Eq. (36) allows to go back step by step obtaining the velocity for each 
possible initial protovorton length. 

As said before, the end of the absorption of fermions must be joined to 
the beginning of the process of stabilization of protovortons. In Fig. 3 the two 
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cases of h = ~ and h — 75 are shown. From these curves the total interval 

5 Id 

of time from the formation of the Kibble string to the birth of the stabilized 
vorton may be obtained according to 

At M =L„f^(£). (38) 

The results are respectively At to t = 3.2 x 10~ 38 sec and 4.5 x lCT 38 sec for 
the two cases of h. 

Therefore, since the time after big bang of the GUT transition is at least 
10~ 36 sec, the delay for the formation of vortons appears to be small so that 
T r ~ ~ Tx and their density Eq.(25) should hold. 

The next step is to analyze what happens to vortons, whose density 
evolves according to Eq.(25), when the universe cools down to the temper- 
ature Tew for the electroweak transition. If the GUT model is such that 
the zero-modes for heavy quarks disappear, most of vortons collapse and the 
non-equilibrium process may allow to produce the needed matter-antimatter 
asymmetry0 . But if all vortons lose these zero-modes instantaneously the 
reheating would be so large that the baryogenesis would result very much 
dilutedcJ . 

The rate of destabilization of vortons depends on the model, and presum- 
ably it is due to the mixture of Higgs mechanisms at GUT and EW scales 
where the latter is responsible for the loss of zero-modes which allowed the 
exotic fermion to be stable inside the string. Therefore it is likely that, being 
the rate of decay of the fermion outside the string ~ olqut rnx, the rate of 
destabilization of vorton is 

/ T \ ^ 

7 ~ Oi GVT m x y-JTTj ~ 1Ql1 secrl ' ( 39 ) 

which is smaller than the Hubble parameter at the beginning of the EW 
transition thus smoothing considerably the reheating effectS . 

Also depending on the model new zero-modes, this time corresponding 
to ordinary fermions, may appear at the EW temperature. Then vortons 
stabilized by them may be formed. But since the loops are collapsing due 
to the loss of the original zero-modes, one has to see whether the ordinary 
fermions succeed in being absorbed to reach the required density before the 
strings disappear. It is understandable that the surviving vortons will be the 
long ones. 
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The string is composed now by superheavy bosons X in its inner core and 
by the EW bosons in the outer part. The length of the inner core is still 
and its width so that the energy per unit length both of the Higgs 
potential and magnetic GUT contributions is ~ m 2 x . The outer ring will have 
a width ~ — !— and consequently a length so that the Higgs potential and 
magnetic weak contributions will be ~ « m 2 x . Regarding the density 
of ordinary fermions to have a new vorton, it must be n £ ~ . 

An argument similar to the one following Eq.(27), and since outside the 
string now n Q ~ T^ w , indicates that 

rii ~ h T£ w At , (40) 

which requires that At is smaller than the collapse time r c . Since for masses 
of ordinary fermions h ~ 10 -3 and r c > m ^ because all the string configu- 
ration including its electroweak part must collapse together to avoid energy 
divergences, Eq.(40) will be satisfied by long strings with iV ~ 10 3 and not 
for the most abundant vortons with N ~ 10. 

According to the discussion of Sec. 2, the ratio of the density of long 
vortons Nl ~ 1000 which allow to produce UHECR and that of short ones 
N ~ 10 Eq.(25) which would give the expected baryogenesis should be 

n L{T) _ 27 28 



n v (T) 



1Q- 27 - 10~ 2S . (41) 



It is a very delicate matter to explain this ratio which cannot be related to 
that of the corresponding Boltzmann factors for vortons but more reasonably 
to the Kibble loops at the beginning of the acquisition of fermions which 
means that their energy is one half of that of vortons of the same length L 
according to Eq.(32), giving a ratio which is still too small. 

But one must remember that, at variance from vortons which are plane, 
the protovortons which originated them may be rough and therefore have 
entropy. Thinking that a protovorton is a chain of objects that can be 
horizontal or vertical, the degeneracy without considering the energy of twists 
is d ~ 2 N , which can enlarge the density of long vortons in such a way that 
Eq.(41) may be satisfied. 

The last important point is to assure that the vorton lifetime for emission 
of an X is at least of the order of the universe age to allow the production of 
UHECR in recent times. The estimation may be done through the bounce 
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instanton where, to simplify things, we consider the vorton as an object in 
the x-y plane so that the Euclidean action Se may be taken as the difference 
of the three-dimensional energy of a tube which contracts emitting an X and 
that of the non-contracted one. The most important contribution to this 
difference will be given by the gradient of the GUT Higgs field <fi so that 

Se~ f(^) 2 6dxdvdz~^PL. (42) 
J \oz J Az 

A<fi ~ nix because it is the change from broken to unbroken vacuum. 5 2 ~ -\- 

m x 

is the section of the inner core of the string where GUT fields are concen- 
trated. Finally Az > m x x , because the length for the variation of the tube 
size must be at least of the order of the Compton length of the emitted 
X particle. Since L ~ Eq.(42) with all the above estimations gives 
Se < Nl and the lifetime for the emission of one X is, being the vorton mass 
m v ~ N L m x , 

T „ — e s z < — *— , (43) 
m v N L m x 

which, as an order of magnitude, ensures r larger than to for N L ~ 10 3 . 

In this way we have shown the feasibility of the mechanism for UHECR 
according to Eq.(19). 

IV. GUT models and possible vortons 

To build a Kibble string it is necessary to break an abelian U{\) symmetry 
different from the electromagnetic U(l). If the model is simply 

U(l) x U{1) - U(l) , (44) 

one such infinite string is stable. But if U(l) is contained in a larger group 
G, a discrete symmetry Z x of the continuous U(l) must remain unbroken 
to avoid that the corresponding monopoles make the string unstable by cut- 
ting it. This depends on the Higgs mechanism and in general the discrete 
symmetry does not survive if the Higgs field corresponds to the fundamental 
representation of G. 

On the other hand the string will become superconducting if a fermion 
acquires mass through a Higgs that winds it. Subsequent phase transitions 
caused by different Higgs fields may produce the disappearance of a previous 
zero-mode. 
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A. SO(10) 

It has 45 generators and one of its maximum subgroups is SU(5) x U(l), 
suitable for string formation. 

There should be four subsequent breakings with the indicated relevant 
Higgs multiplets 

SO (10) 577(5) x [7(1)^577(5) x Z 2 (45) 

-&SU(3) C x SU(2) L x U(1) Y x Z 2 ^SU(3) C x U(l) x Z 2 . 

The decomposition in SU(5) x U(l) of the involved Higgs and fermion 
multiplets is 

45 = 1° + 10 2 + TO" 2 + 24° (46) 
126 = l 5 + 5 1 + 10 3 + I5~ 3 + 45" 1 + 50 1 
10 = 5 1 + 5 -1 
16 = 5^ + 10^5 + i-l . 

The first breaking must be done by a 4> A5 non vanishing in the 1° com- 
ponent to keep the 5/7(5) x U(l) invariance. The second by a 126 = ^ to 
preserve SU(5) symmetry. Regarding fermions, since 16 x 16 = 126+120+10, 
only the second breaking produces mass for one of them, the ur, through a 
Majorana term which is invariant under £7(1) as is seen from Eq.(46). There 
also stable infinite strings!! are produced with winding number n = 1 

015 r^oo ' Vu ' A e ,.^00 > jx - , (47) 

06 T 

and loops will be classically stabilized as vortons by the superconducting 
current of vr. Since SU(5) is still valid, which will presumably require SUSY, 
rj d ~ 10 16 GeV. 

The third phase transition does not affect the string because it is done 
by a 04 5 in 24°, with components which break SU(p) but keep the invariance 
under SU(3) C x SU{2)l x U(1)y, that does not feel the U(l) charge e. 

On the contrary the last breaking, the EW one, has several consequences. 
Since ordinary fermions get Dirac mass from the 577(5) products 

5 x 10 = 5 + 45 , 10 x 10 = 5 + 45 + 50, (48) 
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the 5 1 part of <p 10 will give mass to the quark u and 5 1 to quark d. Due to 
the fact that in SO(10) these masses come from a single mass term 16 x 16, 
to avoid the degeneracy of u and d either one takes different expectation 
values in the neutral components of 5 1 and 5 , or we repeat this term twice 
with ad-hoc coefficients, one multiplied by 1O nontrivial in 5 1 and the other 
with cf) 10 nontrivial in 5 1 with the same expectation value. 

Regarding the string, </> 10 cannot be constant everywhere because the U{1) 
charge would give a contribution to the covariant derivative Dq producing 
a divergent energy. Therefore we must accept the possibility of winding </> 5 i 
and </>g-i and also the inclusion of a neutral Z contribution 



A>05i = ~aE ~ ieA e ~ l 9v Z e </v > ( 49 ) 




1 d 



If we assume that the winding number of 5 i is m, the condition to cancel 
its covariant derivative is 

m - \ - X = , (50) 
5 

where x is the contribution of Z. Obviously, to cancel also the covariant 
derivative of (p^-i , which behaves as 0gi , the same condition Eq.(50) holds 
meaning that it would wind with — m. 

But the minimization of the magnetic energy of Z requires m = so that 
4> 5 i and </>|pi do not wind and no zero-modes appear for ordinary fermions. 

Moreover (f) w gives rise to a coupling V R u L and, since 1O does not wind, 
this small massive contribution to the state of pr makes the corresponding 
zero-mode disappear so that all the vortons would collapse at the EW phase 
transition^ . 

It must be noticed that the situation would change! if an additional Higgs, 
which does not generate mass, 210 = 5~ 4 + 5 + .... with expectation value in 
5~ 4 is present. To compensate the covariant derivative of such a nonwinding 
field 

D e (p 5 -4. = -— + iAeAg - ig^Ze </> 5 - 4 , (51) 



r d9 

Zq needs to behave Z e — >• | j- i which introduced into Eq.(49) requires 
m = 1. In this rather artificial way, ordinary fermions would have zero- 
modes and GUT vortons should not collapse. 
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B. Eg 

It has 78 generators and one maximum subgroup is SO (10) x U(l) with 
the relevant chain of breakings 

E 6 -^SO(10) x U(l)-v-+SO(10).... (52) 
...SU(3) C x SU(2) L x U(1) Y 27 > SU(3) C x ?7(1) . 

The descompositions in 5*0(10) x U(l) of the fundamental and adjoint 
representations are respectively 

27 = l 1 + 1(H + 163 (53) 
78 = 1° + 16~3 +163+45° . 

The first breaking must be done by (p 78 = 0^ to keep the invariance of 

SO(10) x 77(1) and since 27 x 27 = 27 5 + 351 5 + 351 A , no fermion gets 

mass in it. Fermion masses come from Higgs in 27 or 351* (which explains 

the decomposition of 126 in Eq.(46)). Thus the second breaking done by 

<f) 27 = to preserve the invariance of 5*0(10) , gives mass to the exotic 

fermions contained in 10~2 of 27 through a term d)-,iib lib i . These 

to ri r 10 - 2 ^ 10 - 2 

exotic fermions, one quark D of charge — | with 3 colours and an electron E 

with its neutrino N may be the carriers of a vorton based on the string 

fo^mie* > A e ^ 1 -- . (54) 

e r 

Due to the fact that 27 is the fundamental (spinorial) representation, a 
Z 2 symmetry is not preserved and the infinite string is not absolutely stable. 
However it is stable enough if the scale r]jj ~ 10 1& GeV for the breaking of 
U(l) is at least one order of magnitude lower than that of breaking of E$. 

In the last breaking, an expectation value of Higgs in 10"^ gives mass to 
ordinary fermions through a term 6 i ib i ib i . Regarding the influence 

J to ^10^2-^163^163 to to 

of 1Q -i on the string, the difference with the case of 5O(10) is that now the 
U(l) charge — |e is common to both 5 1 and 5 1 of SU(5) x U(l), i.e. 

fid e— 

De ^ 1 = \ rd9 + l 2 Ae ~~ ieA(> ~~ igipZe I ' ^ ^ 



fid e— 

Dgcp^i = I -— + i-A e + ieAg + ig^Zg ] (b^-i . (b) 
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If <fi 5 i winds at the EW scale as <p 5 i — > r] EW e ■ , to avoid energy divergence 
coming from Eq.(55a) it must be 

1 . . 

m + - - X = , (56) 

where now \ is the contribution from Aq and Zq. 

If the minimization of energy favours x = ~ §> f rom Eq.(56) m — — 1 
and since tp s i gives mass to u and v they will have zero-modes and the 
corresponding carriers will run in the direction — z of the string axis. Now 
from Eq.(55b) with X = ~\ the energy divergence is avoided if (j^-i does 
not wind and, since it gives mass to d and electron, they will have no zero- 
modes. Conversely if energy minimization favours x = \i ( t ) -^- 1 wm wind as 
— ► rj EW e imd with m — — 1, d and e will have zero-modes and they will 
be the vorton carriers running along — z of string axis. 

Therefore, with a scheme based on Eq, vortons formed at scale rj-y with 
exotic fermions E and D would remain essentially stable down to scale ri EW 
where would incorporate new carriers either u and v or d and e. 

But if the EW breaking is done with a Higgs in 27 of E e which has not 
only expectation value in its component 10~2 but also in 16*, there may be 
an ordinary-exotic mass term0 6 i ib i ib i , which will mix E and D with 

e and d. Being 10 = 5 1 + 5 , 5 1 contains D, E and iV and 5 D, E and N. 

Therefore if the situation that minimizes energy is the first described 
above, i.e. x = ~ \ > d an d e will have no zero-modes and their mixing with 
D and E will destroy the zero-modes of the latter producing our baryogenesis 
mechanism. On the other hand, the new zero-modes of u and v would be 
responsible for the quasi-stability of the long vortons from the EW age to our 
days giving rise, through quantum decay, to a possible source of UHECR. 

This is in fact what happens. To cancel the large distance covariant 
derivative of <\> x -\ Eq.(55) avoiding the divergence of energy, one needs, 
together with the already existing Aq of Eq.(54), 

A e - % - , Z e - ^ l - , (57) 
e r g v r 

with c\ + c 2 = =f| which would allow to choose |ci| = |c 2 | = j in both cases 
to minimize the magnetic energy | (B 2 + B\ \ . 
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But we now have the additional condition of cancelling the covariant 
derivative of 1 . Because of the decomposition of 16 Eq.(46) in represen- 
tations of SU(5) x (7(1) , the mixing term D d may come from 5 1 x 10~5 
and, due to the fact that 5 x 10 = 5 + 45 , invariance is obtained multiplying 

_3 

by a with expectation value in 5 2 i.e. 0g§'05- 1 V' lo -£- To respect the 
symmetry SU(3) C x U(l) the nonvanishing component of (p_3 must be the 
neutral one so that its covariant derivative is 

( 1 d ~& 3 \ 
= \rdO~ i l' A *~ % 2 ZAe ~ ig '* Zd ^ ' (58) 



Assuming a behaviour cp_3 — > e , the cancellation of energy divergence 

5™ 

requires 

13 

m - - - -ci + c 2 = , (59) 

because = — ^ as ^ corresponds to 5 and g'^ to 5. 

With the condition appropriate for our model c\ + c 2 = — \ , Eq.(59) is 

satisfied with c\ = —0.3 , c 2 = —0.2 and m = producing less magnetic 

energy than that of the alternative case c\ + c 2 = \ which would require 

much more difference between |ci| and |c 2 | . 

— — - _>> 

The other mixing term d D comes from 5 2 x 5 times a nonvanishing 1 2 

of 6 1 with a covariant derivative 

/l d 1 5_~ \ 

where Z does not appear because (f> § has no interactions of SU(5) and 
therefore of the SM. If (f> 5 winds with m , to cancel Eq.(60) 

1 2 

m-i + ^c 1 = . (61) 

and it is easy to see that the case \ — ~~ \ gives again more balanced values 

of Ci and c 2 than for \ — +| , but less than those emerging from Eq.(59). 

Then, since it is not possible to cancel simultaneously Eqs.(58) and (60) 

energy saving suggests that only 0_3 will be nonvanishing. 

52" 

For the model of vortons based on E 6 it is better if the breaking of 50(10) 
in the missing part of the chain of Eq.(52) does not follow Eq.(45) but the 
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alternative@ maximum subgroup SU(A) x SU(2) L x SU(2) R to avoid the 
formation of strings U(l) which would complicate things and also to make 
the GUT unification avoiding SU(5). 

V. Conclusions 

If the explanation of UHECR requires the top-down mechanism, this will 
imply new physics beyond the SM. In the case of the superheavy quasi- 
stable particles, possibly indications on the hidden sector of supersymmetry 
breaking will be given. As for the explanation through cosmic strings GUT 
will be explored and, particularly with vortons, details of the groups and 
Higgs breakings may be revealed. 

We have shown the feasibility of this last mechanism. Obviously, due to 
the semiquantitative evaluation of several effects, only orders of magnitude 
could be given. More accurate analysis must be performed particularly re- 
garding the quantum decay of vortons and precise details of the replacement 
of exotic fermions by ordinary ones at the EW transition as carriers of super- 
conducting current inside the string, as it has been proved to be favourable 
in the case of the Eq symmetry. 

Since the features of the UHECR coming from the two quoted alternatives 
of top-down mechanism are similar due to the fact that they are based on 
quark hadronization, a deep theoretical study of the models will be useful to 
falsify some of them. 

As for the evidence of sources based on superheavy microscopical objects, 
particles or vortons, there are analysisi! to see if there is a possible anisotropy 
of observed events towards the massive concentration in the galaxy. To de- 
termine this, as well as the suggested hard component^ of the spectrum 
beyond the GZK cutoff and the required abundance of primary gammas and 
neutrinos, a much larger statistics is needed as will be supplied by the Auger 
project^ . 

Acknowledgments 

LM thanks CONICET, project PICT 0358, and MO CLAF and Fun- 
dacion Antorchas, project A-13798/1, for partial financial supports. 



21 



References 



X C.T. Norman, D.B. Melrose, A. Achterberg, Ap. J. 454 (1995) 60; A.V. 
Olinto, |astro-ph/0002006| ; N. Nagano, A.A. Watson, Rev. Mod. Phys. 72 



(2000) 689. 

2 K. Greisen, Phys. Rev. Lett. 16 (1966) 798; G. Zatsepin, V. Kuzmin, JETP 
Lett. 4 (1966) 78; F. Aharonian, J. Cronin, Phys. Rev. D 50 (1994) 1892. 



V. Berezinsky, |astro-ph/ 0001 163" ; P. Bhattacharjee, G. Sigl, Physics Rep. 
327 (2000) 109.. 

4 J. Ellis, G. Gelmini, J. Lopez, D.V. Nanopoulos, S. Sarkar, Nucl. Phys. 
B 373 (1992) 339; V.A. Kuzmin, V.A. Rubakov, Yadernaya Fisika 61 
(1998) 1122; V. Berezinsky, M. Kachelriess, A. Vilenkin, Phys. Rev. Lett. 
79 (1997) 4302; K. Hamaguchi, Izawa K.- I., Y. Nomura, T. Yanagida, 
hep-ph /9903207 ; K. Hamaguchi, Y. Nomura, T. Yanagida, Phys. Rev. D 



59 (1999) 063503; R. Benakli, J. Ellis, D.V. Nanopoulos, Phys. Rev. D 
59 (1999) 047301; L. Kofman, A. Linde, A. Starobinsky, Phys. Rev. Lett. 
73 (1994) 3195; D.J.H. Chung, E.W. Kolb, A. Riotto, Phys. Rev. D 59 
(1999) 023501; V.A. Kuzmin, I.I. Tkachev, JETP Letters 69 (1998) 271. 

5 C.T. Hill, D.N. Schramm, T.P Walker, Phys. Rev. D 36 (1987) 1007; P. 
Bhattacharjee, N.C. Rana, Phys. Lett. B 246 (1990) 365; P. Bhattachar- 
jee, C.T. Hill, D. Schramm, Phys. Rev. Lett. 69 (1992) 567; G. Sigl, Space 
Sc. Rev. 75 (1996) 375; G. Sigl, S. Lee, D.N. Schramm, P. Bhattacharjee, 
Science 270 (1995) 1977; R.J. Protheroe, P. Johnson, Astropart. Phys. 4 
(1996) 253; R.J. Protheroe, T. Stanev, Phys. Rev. Lett. 77 (1996) 3708. 

6 T.W.B. Kibble, J. Phys. A 9 (1976) 1387. 

7 P. Bhattacharjee, Phys. Rev. Lett. 81 (1998) 260. 

8 V. Berezinsky, A. Vilenkin, Phys. Rev. Lett. 79 (1997) 5202; V. Berezinsky, 
P. Blasi, A. Vilenkin, Phys. Rev. D 58 (1998) 103515. 

9 E. Witten, Nucl. Phys. B 249 (1985) 557. 

10 L. Masperi, G. Silva, Astropart. Phys. 8 (1998) 173; L. Masperi, M. Orsaria, 
Nucl. Phys. B (Proc. Suppl.) 75 A (1999) 362. 



22 



n L. Masperi, M. Orsaria, Int. J. Mod. Phys. A 14 (1999) 3581. 
12 M. Takeda et al., Phys. Rev. Lett. 81 (1998) 1163. 

13 G. Sigl, S. Lee, D. Schramm, P. Coppi, Phys. Lett. B 392 (1997) 129 
Z. Fodor, S.D. Katz, |astro-ph/ 0008204] ; L. Masperi, M. Orsaria, |astro- 



ph/0008525|. 



14 Y. Chikashige, J.i. Kamoshita, |astro-ph/ 98 12483] , 

15 Yu.L. Dokshizer, V.H. Khoze, A.H. Mueller, S.Y. Troyan, Basics of Pertur- 
bative QCD (Editions Frontiers, 1991). 

16 M. Birkel, S. Sarkar, Astropart. Phys. 9 (1998) 173. 

17 M. Blanton, P. Blasi, A. Olinto, |astro-ph/ 0009466. 



18 O.E. Kalashev, VA. Kuzmin, D.V. Semikov, |astro-ph/ 0006349] ; S.L. 
Dubovsky, PG. Tinyakov, JETP Lett. 68 (1998) 271. 

19 G. Medina Tanco, |astro-ph/0009336 . 



20 B. Carter, |hep-ph/99095T3] . 

21 B. Carter, A.-C. Davis, hep-ph / 9910560. 

22 R. Brandenberger, B. Carter, A.-C. Davis, M. Trodden, Phys. Rev. D 54 
(1996) 6059. 

23 C. Martins, E. Shellard, Phys. Rev. D 57 (1998) 173. 

24 W.B. Perkins, A.-C. Davis, Phys. Lett. B 393 (1997) 46. 

25 A.-C. Davis, S.C. Davis, Phys. Rev. D 55 (1997) 1879. 

26 S.C. Davis, A.-C. Davis, W.B. Perkins, Phys. Lett. B 408 (1997) 81. 

27 J.L. Rosner, [hep-ph/9907438] . 

28 P. Langacker, Phys. Rep. 72 (1985) 125. 

29 V. Berezinsky, AA. Mikhailov, |astro-ph/9810277| ; GA. Medina Tanco, 
A. A. Watson, |astro-ph/9903T82j 



23 



A.M. Hillas, Nucl. Phys. B (Proc. Suppl.) 75 A (1999) 109. 

Pierre Auger project design report (1997), [http: / / www-td-auger. fnal.gov: 



24 



26.0 



25.5 - 




"1 1 1 1 I 1 1 1 1 

18.5 19.0 19.5 20.0 20.5 

Log ( E [eV] ) 



Fig.l. Experimental UHECR including • AGASA and A Fly's Eye data together 

- 3 2 

with a fit corresponding to a galactic flux ~ E followed after the ankle by a possible 
extragalactic flux ~ E " 2 8 which exhausts beyond the GZK cutoff giving way to the hard 
component of vortons in halo ~ E" . 
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Fig.2. Fit of calculated energy spectrum of UHECR from vortons uniformly distributed 

in universe ~ E " <1+K) without considering interaction with CBR resulting in a redshift 
correction K = 0.01 1 . 




Fig.3. Velocity of loop contraction v before reaching vorton size L v . Continuous line 
represents protovorton stabilization. Dashed and dash-dotted lines give protovorton 
formation with carrier absorption coefficient h = 1/5 and 1/13 respectively. 



